Scattering of electrons from an interacting region 
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We address the problem of transmission of electrons between two noninteracting leads through 
a region where they interact (quantum dot). We use a model of spinless electrons hopping on a 
one-dimensional lattice and with an interaction on a single bond. We show that all the two-particle 
scattering states can be found exactly. Comparisons are made with numerical results on the time 
evolution of a two-particle wave packet and several interesting features are found. For N particles 
the scattering state is obtained within a two-particle scattering approximation. For a dot connected 
to Fermi seas at different chemical potentials, we find an expression for the change in the Landauer 
current resulting from the interactions on the dot. We end with some comments on the case of 
spin-1/2 electrons. 
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An understanding of the behavior of electrons inter- 
acting with each other in a localized region has been a 
challenging problem in theoretical physics. Recently it 
has attracted much attention in view of the experimental 
interest in transport across quantum dots and the Kondo 
effect in a quantum dot As a prototypical model, let 
us consider two ideal leads, where all electronic interac- 
tions can be neglected, connected to a region (a quantum 
dot) where the electrons interact. One is interested in the 
current through the dot in response to an applied voltage 
difference between the leads. 

As has been discussed in Ref. Q, there are several 
different but equivalent theoretical approaches. In the 
nonequilibrium Green's function (NEGF) approach the 
initial density matrix, of the two reservoirs (taken as ideal 
Fermi liquids in equilibrium at different chemical poten- 
tials) and the dot (in an arbitrary initial state), is evolved 
in time. The coupling between the reservoirs and the dot 
is switched on adiabatically and one looks at steady state 
properties of the resulting density matrix. A related ap- 
proach is the quantum Langevin method where the reser- 
voirs are treated as sources of noise and dissipation. A 
second approach is to view this as a time-independent 
scattering problem and to look for many-particle scat- 
tering states which have the correct asymptotic form in 
the leads. This is in the spirit of the Landauer formalism. 
In the case where there are no interactions in the dot re- 
gion, exact results for the current and other steady state 
properties can be obtained, and all three approaches give 
identical answers [H, 0, 0, 0] • 

The interacting case however is much more difficult to 
study. For a single dot connected to noninteracting leads, 
some results using the NEGF method have been obtained 
using the so-called non-crossing approximation For 
an integrable model, namely the interacting resonance 
level model, Mehta and Andrei used the scattering ap- 
proach to solve the problem exactly . Using the Bethe 
ansatz, they were able to express all A^-particle scatter- 
ing states in terms of the two-particle S-matrix which 



is known exactly. They considered a continuum model 
with a linear spectrum which makes it integrable. The 
iV-particle scattering matrix for electrons interacting in 
a quantum dot has also been studied in Ref. @, 0] • 

In this Letter, we study a lattice version of the model 
considered in Ref. Q. We show here that using the 
Lippman-Schwinger method all two-particle eigenstates 
of this model can be found exactly. The form of the S- 
matrix indicates that the model is not solvable by the 
Bethe ansatz. We examine the S-matrix and compare 
it with numerical experiments on scattering of a two- 
particle wave packet. We also study many-body trans- 
port in this system by considering TV-particle states cor- 
responding to left and right leads with different chemical 
potentials. We obtain an expression for the change in the 
Landauer current arising from the interactions. 

We note that the study of two-particle scattering states 
is in itself of interest [9|, llfj , apart from being the starting 
point for the study of many-particle states necessary to 
understand transport. Recently, Goorden and Biittiker 
ll| have studied a set-up with two disconnected conduct- 
ing wires and with electrons in the two wires interacting 
weakly in a localized region. Using first order perturba- 
tion theory, the two-particle S'-matrix was evaluated and 
used to extract information on transmission and corre- 
lations in a two-particle scattering experiment. In our 
single channel case, we will show that the antisymme- 
try of the wave functions leads to striking asymmetries 
in the S'-matrix. In another interesting recent work, the 
S'-matrix in a model of two photons interacting with a 
localized atom was studied 



We consider a tight-binding one-dimensional lattice 
with spinless electrons. The model considered describes 
an interacting dot on the sites x = 0, 1 which is con- 
nected to two noninteracting one-dimensional leads on 
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either side. The Hamiltonian is given by 

H = H L + H D + V c , where 



(1) 



H, 



X— — oo 



(4^+1+4+1^)' 



= -7( c o c i + 4 c o) + e (™o + «i) + Unom, 

and Vc = -7 , (cL 1 c + 4c_i + c|c 2 + 4 c i): 

where n x = 4 c z is the number operator at site a;, and 
implies omission of a; = —1, 0, 1 from the summation. We 
set the lattice spacing and h to 1. In this paper we only 
consider the case 7 = 7' = ! and e = corresponding 
(for U = 0) to the case of a perfectly transmitting dot 
but the general case can be treated similarly [l^ . 

Scattering states: We first show how one can obtain 
all the two-particle energy eigenstates exactly for this 
problem. Consider the noninteracting Hamiltonian Hq = 
H with U = 0. For this case, the one-particle eigenstates 
have the form <pk{x) = e zkx with energy E% = — 2cosfc, 
where —ir<k<Tr. Now consider a two-particle incom- 



ing state given by 0k (x) = e 



i(k\xi+k 2 x 2 ) p i(k 2 xi-\-kix 2 ) 

with k = (fci,fc2) and x = {x\,X2)- The energy of this 
state is E^ — E^ + Ef. 2 . A scattering eigenstate of 
iJ = Hq + V (where V = Urioni) with energy E is related 
to a state |0) of Hq by the Lippman-Schwinger equation 

|V) = |0) + G+{E)V\i>), (2) 

where Gq^-E 1 ) = 1/(E — Hq + ie). In the two-particle 
sector, in the position basis |x) and with an incident state 
< x|0) = k (x), Eq. © gives 



^k(x) = k (x) + UK Ek (x) Vk(0), 



(3) 



where J^Jx) = (x|G^(E k )|0) and = (1,0). 
We can determine V'k(O) using Eq. ([3]), -0k (0) = 
0k(O)/[l — UKE k (0)]. The matrix element Ke^(x) is ex- 
plicitly given by 

K E k (x) = gt k (xi - l,x 2 ) - g%Jxi,X2-l), (4) 

where = [1/(2tt) 2 ] d^d^e*** / ( E k - 

Eq+ie) is the usual two-dimensional lattice Green's func- 
tion. It is instructive to look at the asymptotic form of 
the scattered wave function [14| ; this can be obtained by 
the saddle point method, the contribution to the integral 
in Eq. coming from the region near D q = E^. Apart 
from a factor Uipk{0), we find asymptotically that 



nt(x) 



(±1 _ j\ e i(k' 1 x 1 +k' 2 x 2 ) 

4tt 1 /2 (r/ro) 1 / 2 



with 



•?'2 



sin(fcj) sin(fc 2 ) 
£J k = - 2cos(^ 
r 



(e - e lfc2 ), 
(5) 

— , where Xij sin(fc-) > 0, (6) 
2cos(/c 2 ), 



T 1 + x 2 



and ro 



2\l/2 

2) j 

sm 2 (A;i)+sin 2 (fe^)] 1 /2 



(7) 

(8) 



where the ± sign in Eq. ([5]) corresponds to E^ ^ 0. The 
antisymmetry of the wave function is implicitly hidden 
in the x-dependence of k'. [The expression in Eq. ([5|) 
is clearly more complicated than the Bethe ansatz would 
have given which is a superposition of only four pairs of 
momenta, namely, (±&i, ifo).] The physical interpreta- 
tion of the above solution is as follows. Two electrons 
with initial momenta (ki,k 2 ) emerge, after scattering, 
with momenta (k'i,k 2 )- Energy is conserved as implied 
by Eq. ([7|). (Momentum is not conserved because the 
interaction term Un^ui breaks translation invariance). 
The velocities of the electrons are given by v\ = 2 sin(fc^) 
and V2 = 2sin(fc 2 ); Eq. ((6]) expresses the fact that the 
electrons observed at (2:1,0:2) must reach there at the 
same time after collision. Note that we can equivalently 
think of this problem as that of a single electron in a 
two-dimensional (2D) lattice moving in the half-space 
x\ > X2, with a hard wall along the diagonal x\ — x 2 and 
a single impurity at the site 0. The particle flux J ■ dS 
in a given direction tan(Q) = x 2 jx\ in the ID problem 
corresponds, in the II? problem, to the rate at which two 
particles are scattered with velocity ratio v 2 jv\ = tan(9). 
Instead of the usual scattering cross-section, it is useful 
here to calculate the scattering rate for unit two-particle 
density at the site 0. This is given by 



\f{0)\ 2 d6 = 
~ I 1 



J-dS 



1 



|0 k (O)| 2 \1/U - K Ek (0)\ 2 
cos(fcJ - k' 2 )\ [sin 2 (fci) + sin 2 (fc 2 )] 



2?r| sin 2 (fci) cos(fc 2 ) + sin 2 (fc 2 ) cos(fc£)| 



d9, 



(9) 



I sin 2 (fci) cos(fc 2 ) + sin 2 (£; 2 ) cos(fc()| ' 



where k[, k' 2 are known in terms of 9. Experimentally it 
may be simpler to find the number of particles scattering 
within an energy interval dEy (energy conservation 
implies that dEy + dE k > = 0). Defining P(Ek 1 ,Ek 2 — > 
E K ,E K )dE K = |0 k 2 (O)| 2 |/(0)| 2 ^, we find that 
P = [l-cos(fci-fo)][l-cos(fc£-/4)]/{|l/Z7 - 
K Ek (0)\ 2 27r|sm(fc' 1 )sin(fc 2 )|}. 

For the two-particle case it is more useful to study wave 
packets. We now consider the time evolution of wave 
packets and see how well the predictions of the scattering 
theory hold. The scattering states given by Eq. §5§ are 
the full set of allowed two-particle energy eigenstates (for 
U > 7r one gets an additional bound state [HI]). These 
can be generated by a unitary time evolution of the un- 
perturbed states which form a complete set. Hence these 
states also form a complete set, and any two-particle wave 
function can be expanded using this basis. Thus the time 
evolution of an initial wave packet ^(x, t = 0) is given 
by 

*( X '*) = (27r)2 J dqi J rf ^2 a(q)^ q (x) e- 1 ^ 4 , 
where a(q) = ^ *(x,i = 0) Vq(x). (10) 
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FIG. 1: Plot of the evolution of an incident wave packet (a) 
after passing through the origin with U = in (b) and (7 = 2 
in (c). Note the strong scattering at an angle 9 — — 7r/4. 




FIG. 2: Plot of the number of particles scattered into a given 
direction for incident wave packets with different energies and 
U — 2. The bold lines show the results from scattering theory 
estimated using |/(#)| 2 and the incident particle density at the 
origin (inset). Inset shows |^I'i nc (0, t)| 2 . 



The time evolution can be studied quite accurately be- 
cause of our knowledge of the exact basis states. In evalu- 
ating the basis states, for small (xi, x 2 ) < 15, we evaluate 
the Green's functions g J (x) exactly using recursion re- 
lations [H| relating these to g^ k (0,0) and (1, 1). For 
larger (xi,x 2 ) we use the asymptotic forms which are 
quite accurate. We find that in our computations the 
normalization of the wave function is preserved to within 
0.5%. In Fig. [T]we show the typical time-evolution of a 
wave packet with initial position and momentum local- 
ized at x = (-3.922, -5.0) and q = (2.554, 0.785) respec- 
tively and with widths 5x « 5q « 1 and E 0.25. These 
initial conditions have been chosen so that the two parti- 
cles reach the site at roughly the same time; this max- 



imizes their interaction. The initial wave packet shown 
in Fig. 1 (a) evolves at time t = 20 to (b) for U = 
and to (c) for U = 2. For the scattered wave function 
in Fig. [1] (c) we can count the number of particles scat- 
tered into a given direction. This is plotted in Fig. [2] for 
incident wave packets with different energies. We also 
compare this with the scattering theory prediction by 
plotting |/(6') | 2 multiplied by the time-integrated inci- 
dent two-particle density at the origin. The comparison 
can be seen to be very good. 

Transport calculation: We will now turn our at- 
tention to quantities of interest in transport calculations. 
The current density is given by the expectation value of 
the operator j x = —i(c\.c x +i—h.c.) in the scattering state 
iV-'k) = |0k) + |<Sk). The current in the incident state is 
given by {4>k\j x \4>k) = 2[sin(fci) + s\n(k 2 )]N , where N is 
the total number of sites in the entire system. The change 
in current due to scattering, Sj(ki,k 2 ) = (V'kl jxlV'k) — 
('/'kljzl'/'k), gets contributions from two parts, namely, 
is = (SkljxlSk) and j c = (Sk|jx|</>k) + {^k\jx\Sk), and 
is of order 1, i.e., it is a factor of N smaller than the 
current in the incident state tldl. We find that 



k 2 ) 



2\M0)\ 2 Im[K E ^0)} 
\1/U - K Ek (0W 



[sgn(ki) + sgn(k 2 )}, 
(11) 

where sgn(k) = \k\/k. 

N -particle scattering state and change in the Landauer 
current: We now consider the problem of calculating 
the current in a situation where the interacting region 
is connected to left and right leads which are at zero 
temperature and chemical potentials \il and fin respec- 
tively. In that case we have to consider an initial state 
with Nl electrons in positive momentum states filling 
1-particle energy levels up to [Xl and Nr electrons in 
negative momentum states filling levels up to fj,R. Let 
7Y = N E + Nr and let us denote this TV-particle inci- 
dent wave by |<j&W) = |kW), where kW = {kik 2 ...k N }. 
One then needs to find the corresponding scattering state 
and compute the particle current. An exact solution 
for the TV-particle scattering state looks difficult. Wc 
will therefore restrict ourselves to an approximation in 
which only two-particle scattering is taken into account. 
Within this approximation, the scattered wave is given 
by IV'kjv) = l ( / , k JV ) + |5'k J v); where the transition amplitude 
to a wave vector q^v = {qiq 2 ---qN} is given by 



(qjvl^k 



- EM 

q2k 2 



with (qalSika) = 



,p+p' 



Q2 



:q2|^k 2 )(q w - 2 |k' / 



(O)0k 2 (O) 



(l/U- 



N-2h 



fl2l 



i^ k2 (0))(£ k2 -£ q2+ie ; 



Here q2 (k^) denotes a pair of momenta chosen from the 
set qjv (kjv), and q' (k') denotes the remaining N — 2 
momenta. P (P') are the appropriate number of per- 
mutations. Using Eq. (|12p. we can calculate the cur- 
rent expectation value for the state IV'kjv) to order U 2 . 
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(At order U 2 , there are also contributions to the current 
from three- and four-particle scattering, but we will ig- 
nore those here). The current in the incident state |</>k N ) 
is given by {<t>\j x \<t>) = 2Ef =1 sinfo)]^- 1 . The cor- 
rect normalization is obtained by dividing by a factor 
J\f N which then gives in the continuum limit: j inc = 
[f fe£ dk 2sin(fc)- /*« dk 2ain(fc)]/(27r) - (fi L -fi R )/{^), 
where kL,R — cos^ 1 (— Hl,r/2), and we have used dk = 
dE/\dE/dk\ = dE/\2 sin(fc)|. Inserting factors of h and e, 
this gives the expected Landauer current / = (e/h)(/j,L — 
Pr) and Landauer conductance G = e 2 /h. The change 
in the Landauer current due to two-particle scattering is 
given by a sum of two-particle currents from all possi- 
ble momentum pairs: Sjx = (1/2) J2 r s $3{kri k s ) N N ~ 2 
which, with the same normalization as used earlier, gives 



Sj 



1 



A' 



2(2tt) 2 



dkidk 2 Sj(ki,k 2 ), 



(13) 



where the integrations are over the full range of allowed 
momenta [— ir, n], and Sj(ki,k 2 ) is given by Eq. (fTTj) [ex- 
panded to order U 2 ]. Using the fact that 5j(k\, k 2 ) van- 
ishes whenever k\ , k 2 have opposite signs and converting 
Eq. (fT3"|) to energy integrals, we find the following cor- 
rection to the Landauer current, 



dE kl / ,1E, 

-2 J flR 



k 2 



- / dE kl / dE k , 



x p(E kl )p(E k2 ) UH\fa lM (0)\ 2 Im[K EhiM (0)], 

(14) 



where p(E) = l/(27iV4 — E 2 ) is the density of 
states. The quantity in Eq. (fl"4|) is negative because 
Im[KE k {0)] < for all values of k. In the zero bias limit 
Mi ~~ * Mfli Eq. ([14]) vanishes as U 2 (pr, — M-r) due to the 
contribution coming from the first set of integrals; thus 
G is less than e 2 /h by a term of order U 2 . 

Finally, let us briefly discuss the case of spin- 
1/2 electrons. We consider the Hamiltonian H = 

-E^-ooE ct =t,i (4 l0 - c *+i,<r + h.c.) + Un^n oi . The 
interaction at the site can cause scattering between 
two electrons in the singlet channel but not in the 
triplet channel. The scattering of two electrons in 
the singlet channel can be studied exactly using the 
Lippman-Schwingcr formalism just as in Eqs. (OOj 
except that the wave function for the state |(/>k) = 

\ki,V,k 2 ,[) = -|fc 2 ,|;fci,T) is now s iven by <Mx) = 

e i(k 1 x 1 +k 2 x 2 ) ^ an( j foe Green's function is given by 
JG^x) = [1/(2tt) 2 ] /* w d qi dq 2 e^/( E* E q + ie). 
Finally, we can argue as in the spinless case, that in the 
presence of a Fermi sea, the scattering reduces the Lan- 
dauer conductance by a term of order U . 

Discussion: We have shown how the Lippman- 
Schwinger formalism can be used to obtain exact results 
for two particles scattering from an interacting region. 
This method can be applied to other cases, such as the 



two- wire system studied in Refs. the case of spin-1/2 
electrons as mentioned above and the case with interac- 
tions on more than one bond. We have demonstrated 
how scattering theory can be used to understand nu- 
merical results for a two-particle wave packet moving 
through the interacting region. Finally, we have con- 
sidered the problem of many-particle transport across 
the interacting region; we find that two-particle scat- 
tering reduces the zero-temperature Landauer conduc- 
tance by a term of order U 2 . This calculation is nontriv- 
ial since it considers many-particle states and is a fully 
nonequilibrium treatment. We expect the two-particle 
scattering approximation to be valid at low densities 

r/tt << 1 [3 since the 3-particle correction given 
by J J J dkidk 2 dk 3 5j(ki,k 2 ,k3) would be smaller by a 
factor of order kh,Rl^- In this paper we have considered 
the simplest case with interactions on a single bond and 
no impurities. For interactions on more than one bond, 
the form of the two-particle S'-matrix would change but 
the qualitative conclusions remain the same [HI]. In the 
presence of impurities however, a term of 0{U) appears 
in the correction to G and this could lead to an enhance- 
ment of G, depending on the sign of U . More gener- 
ally, interactions can lead to dephasing and suppression 
of weak localization thereby increasing G [13] ■ 
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